Magneto-Roton Modes in the Ultra Quantum Crystal: Life and Death at Large 

Magnetic Fields 
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The Ultra Quantum Crystal phases observed in quasi-one-dimensional conductors of the Bechgaard 
salts family under magnetic field exhibit both Spin Density Wave order and a Quantized Ifall Effect, 
which may exhibit sign reversals. We examine the case with no sign reversals. As the field increases, 
Magneto-Roton modes , which have minima at quantized wave vectors, merge, or not, in the single 
particle excitation continuum. In the high field phase, the Hall conductivity at zero temperature 
is zero; at the transition from the A'^ = 1 phase, it has two Magneto-Roton modes, the energies of 
which increase rapidly with magnetic field. The modes survive until a critical field above which the 
high field phase looks more like a perfect nesting insulating Spin Density Wave. 
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I. INTRODUCTION 

. Organic conductors of the Bechgaard salts family, 
{TMTSF)2X where TMTSF = tetramethylselenafulva- 
lene are quasi-one-dimensional (quasi-lD) systems, which 
have been found over the last few years to exhibit fasci- 
nating properties under magnetic field [0-|5j. The typi- 
cal hierarchy of their transfer integrals is: ta — SOOOiir, 
tb = 300i^, tc = \QK . In three members of this family 
{X — CIO4, PFq, i?e04), the metallic phase is destroyed 
by a moderate magnetic field H applied along the c direc- 
tion, perpendicular to the most conducting planes (a, 5). 
A cascade of magnetic phases, separated by first order 
transitions appears as the field intensity is stepped up: 
within each sub-phase, which are called "Ultra Quantum 
Crystal" (UQC in the following) [0, a Field Induced 
Spin-Density Wave Phase (FISDW) is stabilized with a 
peculiar electronic structure, characterized by a small 
number of exactly filled Landau levels (bands in fact) [Q . 
The Landau bands are separated by a hierarchy of gaps 
6n which oscillate with the magnetic field The phase 
labelled by N has a Hall conductivity a^y — 2Ne'^/h, and 
is characterized by the largest gap (5 at at the the Fermi 
level. A sum rule = A^, connects all the gaps to the 
order parameter A. Thus each UQC sub-phase exhibits 
a Quantized Hall conductivity, which is the first example 
of a Quantum Hall Effect in a 3D system. This cascade 
of quantized phases result from an interplay between the 
nesting properties of the Fermi Surface (FS), and the 
quantization of electronic orbits in the field: the wave 
vector of the FISDW varies with field so that unpaired 
carriers in a subphase are always organized in completed 
filled Landau bands. As a result, the number of carriers 
in each subphase is quantized, and so is the Hall conduc- 
tivity 0,^. The condensation of the UQC phases results 
from the peculiar electronic structure of open Fermi Sur- 
face metal under magnetic field: because of the Lorentz 



force, the electronic motion becomes periodic and con- 
fined along the high conductivity direction of the chains 
(a direction) ||^ . The periodic motion of the electrons in 
real space is characterized by a wave vector G = eHb/h, 
and an orbital energy huic — vpG/2, b being the inter- 
chain distance. (In the rest of this Letter, wave vectors 
will be expressed in units of G). 

As a result, the static bare susceptibility of the normal 
phase, Xo(Q) can be expressed as a sum over weighted 
strictly ID bare susceptibilities which diverge at quan- 
tized values of the longitudinal component of the wave 
vector Q" = 2kp+n The largest divergence signals 

the appearance of a SDW phase with quantized vector 
Q|| = 2kF + N. This Quantized Nesting Model (QNM) 
describes most of the features of the phase diagram 
in a magnetic field. It has been shown recently to ex- 
plain the experimental observation of the Hall plateaux 
sign reversal when the field varies [||-||. Most plateaux 
exhibit the same sign. (By convention we will refer to 
these plateaux as positive ones). 

It has been pointed out some years ago [|l^,0 that 
FISDW were a novel kind of electron-hole condensate, 
an Ultra Quantum Crystal: their condensation energy 
is smaller than the orbital energy huc = hvFG/2. As 
such, it exhibits specific collective modes with a Magneto- 
Roton (hereafter MR) structure which has its minimum 
at a value of the parallel component {q\\) of its momentum 
equal to nG {n integer). The original paper on the UQC 
collective modes suggested that MR modes might exist 
at (7|| = 2G. 3G, etc., but only the mode at G was proved 
to exist [00]. 

Recently one of us (PL) remarked that the UQC with 
"negative " Quantum Hall Number, the "Ribault Phase" , 
exhibits a MR spectrum with no analog in Condensed 
Matter Physics, different from that of the majority sign 
UQC phases iQ. Indeed, when the sequence of observed 
Quantum HaU Numbers [§ is : TV = 1, 2, -2, 3, 4, 5, 6, 7, 
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the Ribault Phase (i. e. N = —2) has MR modes with 
g|| = 4 or 5 well within the single particle gap, and no low 
energy mode at gn = 1 or 2. It was shown in that the 
Ribault phase "contaminates" the neighbouring phases 
in parameter space, as well as those thermodynamically 
close by. In the latter example, the phase N — 2 (resp 
N ~ 3) has at least one additional mode at g|| — A (resp. 
(7|| = 5). Examples of Hall sequences with more sign 
changes have been shown to have more structure yet in 
the collective modes, with possibly 4, etc., MR modes 
within the single particle gap. Furthermore, in contrast 
to the "usual" MR, the two MR in the Ribault phase 
have large monotonic variations, with opposite signs, as 
the field varies (l2j . 

The latter findings prompted us to reexamine in more 
details, and in a quantitative way, the MR of the usual 
UQC, when no sign change of the Hall Effect occurs as 
the field varies jo) . 

This Letter gives for the first time numerical esti- 
mates, based on the RPA, of the MR modes energies in 
(TMTSF)2X , of their magnetic field dependence, and 
of the MR transverse effective mass. We find a distinctly 
different behaviour of the UQC at high fields {N = 1 or 
0) from that at low field > 2. 



II. MR IN THE GENERIC CASE 

We are interested in the usual UQC with no sign inver- 
sion of the Quantum Hall Effect. Therefore it is sufficient 
to consider the following electronic dispersion relation: 

eik)^VF{\k,\-kF) + e^{k^), (2.1) 
e±{kj_) = —2tb cos kyb — 2tc cos kzC — 2t^ cos 2kyb 

e(fcj^) is a periodic function which describes a warped FS. 

The normal metal-FISDW instability hnc Tcn{H) is 
given by: 

Xo(Q,T,^,i/) = ^nll{Q±)xV'{Q\\-n,T,N) = 1/A 



(2.2) 



A is the electronic interaction constant. Eq.(2.2) ex- 
hibits the structure of xo a-s the sum of one dimensional 
terms Xo^ shifted by the magnetic field wave vector 
G = eHb/h. xl" « -ln(max{t;j.(2fcj. - q),T)/eF). In 



(i. e. when = 0). Now define a generalized instability 
temperature TN±m{^q\_)'- 
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In (2.4), we have used the sum rule: = 1. For m — 



ana q± = 0, TN±m — Tcn, the ordering temperature for 
the A^th subphase. For m ^ 0, TN±m{q±) generalizes 
the definition of the critical temperatures on either side 
of phase A^ in the {T,H) plane. TN±m{q±) are at most 
equal to the virtual transition lines Tf^±m which can be 
drawn in the A^th subphase part of the phase diagram 
and which represent virtual transition lines to phases 
with slightly larger free energy than the A^th subphase 
I p^ . In the (T, H) plane, there is an infinite number of 
continuous lines crossing the phase diagram. The upper 
limit of this family is the actual (continuous non ana- 
lytic) transition line from the normal metal to the UQC; 
this line coincides piecewise with the transition lines la- 
belled by the successive integers describing the Quantum 
Hall conductivity. See Fig.(l) . An example of computed 
network of transition lines was given in pq | 

MR energies have been derived within the RPA ||l^,|ll| , 
by looking at the poles of the spin-spin correlation func- 
tion of the ordered phase. The electronic orbital motion, 
together with the RPA treatments of interactions, results 
in effective scattering potential energy terms which cou- 
ple electron states with wave vector k\ | and k\\+ 2kF + n, 
(n integer) on either side of the Fermi Surface. This 
results in a series of gaps in the condensed phase disper- 
sion relation ||^, corresponding to the various potential 
scattering terms. The simplest approximation which cap- 
tures the esssential physics resums to all orders the gap 
<^7V = ^Af X A at the Fermi level and takes all other gaps 
into account to second order in perturbation [ p^|jTl| ] . The 
limit 2bnl^c <C 1 is assumed to hold; in fact this ratio 
is small at < 4T, and is small again at large 77, since 
28* = 3.53 X T* is a fixed bound. Then the equation for 
collective modes in the UQC phase A^ reduces to: 



In 



eq.(2.2), the coefficient /„ depends on the dispersion re- 
lation and H: 

IniQ±) = expi[{T±{p + Q±/2) +T±{p- Q±/2) + np] (2.3) 
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is defined in eq.( p.4D , and g|| = 2fci? -t- 
<< 1 . Xo and Fo are for n = the 



where T±{p) = (l/huic) /q £±{p')dp' and A{p) denotes 
the average of A over p. 

Define T* as the FISD W o rdering temperature 
in infinite field. From eq.(2.2) , we have T* = 
(27/7r)i?o exp(— 1/A). (7 is Euler's constant and Eq a 
high energy cut-off). T* is equal to the ordering temper- 
ature for the perfectly nested Fermi Surface at any field 



where T/vir 

m + S, with 6 
objects discussed in J]l6| in connections with collective 
modes of SDW. Eq.(|^) finally yields, setting x^.m — 
w(g|| =m + S, qj_)/l2SN), with {q^ tt/c): 
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We have computed 



.,1 J / \ roo 7 tanh(-:7^ cosh n) 

With ft x) = „ du W 2-^. 

^ ^ jO cosh^ u—x^ 

MR energies at OK i n the generic UQC phase labeUed by 
N = A from eqs.( |2.4| ), varying q± in order to get the locus 
of the various MR minima in the q± direction. As a by- 
product, we obtain for the first time a numerical estimate 
of the transverse effective mass m± defined in phase N 
by u;N{q±,q\\ = m) = w]^™ + q\/2m± {fi = 1). The pa- 
rameters we took for the numerical estimates are relevant 
to the Bechgaard salts, and have been used previously in 
the literature t„ = SOOO/sT, h = 265K, t'^ = lOK, 
T* = 6.44. (Those parameters may change with pres- 
sure, notably t'f^). The results are shown on fig. (2). We 
find two modes, at = 1 or 2 which vary little with H, 
and a third mode, rather close to the upper edge of the 
gap, at q\i = 6,q± — 0. This mode is yet another ex- 
ample of "contamination" by the Ribault phase N = —2 
discussed in ref. ||l^ . ft has a monotonic H negative vari- 
ation, and drowns into the single particle continuum for 
H < 3.75T. We have looked for MR mode with g|| =8, 
in phase 4, which might arise from contamination of the 
near-by Ribault phase TV = — 4, and found it did not ex- 
ist in the absence of additional terms in eq. ( p.lD , but this 
mode will certainly become alive under pressure [0,0 . 

The RPA evaluation yields a MR minimum ~ .4, 
located at gj. = .0178/6 for H = 3.6T. We find' a very 
large mass anisotropy m±/m\\\N^5^m=i — -9 x lO-^, with 
TO|| = hujmin/vp. The fact that x^^i < 1/2 makes this 
MR mode a good candidate for the kind of specific heat 
anomalies discussed in The mass anisotropy is 

larger for the mode g|| = 2, we find ~ 2.25 x 10^. See 
fig. (3). Note that the MR energy w]^™ is discontinuous 
at the transition from iV to ± 1. This contributes (a 
small term, which corrects the mean field estimate) to 
the transition latent heat On the other hand, w™™ 
goes continuously at the transition to N 
(i. e. a distinct MR with the same energy!) 
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III. LIFE AND DEATH AT LARGE MAGNETIC 
FIELD 

The situation changes qualitatively as soon as the mag- 
netic field exceeds 6.1 T, the critical value for the transi- 
tion from UQC phases 2 to 1. The results are shown on 
fig. (4). The two MR which exist in the low field part of 
phase iV = 1 have different field dependences. The high 
energy mode, with q^ = 2 increases its energy monotoni- 
cally with field, at a rate ~ .6K/T, then merges into the 
single particle continuum for H ^ 7 T. Then only one 
MR at g|| =1 survives, at lower energy, and with a non 
monotonic field dependence. In fact the rate of variation 
of wj^i™ is also large, ~ -1.8K/T for H > 7.5 T. For 
H < 6.5 T, duj'^'J^/dH ~ AK/T. 

At the transition from = 1 to 0, at _ff ~ 8.25 T, two 
MR modes appear discontinuously. The MR at g| | = 1 is 
21% higher in energy in phase 0, and that at gn = 2 is 
close to the single particle continuum, at .95 x 25^. In fact 



this latter mode is also connected to the "contamination" 
of the Ribault pha se a t N = —2, since it connects T2 and 
T_2 through eq. (2.4) [Q. Both MR energies increase 
rapidly with H, at a rate ~ .5K/T for the low energy 
mode, .8K/T for the high energy one. The last surviving 
MR dissolves into the single particle continuum for H > 
12.5 T. For larger H the only collective modes left are 
the usual Goldstone modes at q\\ <C 1. Then the only 
specific features of the UQC left are small single particle 
energy gaps, at most an order of magnitude smaller than 
the FS gap, at large energies hujc, 2h,uj(., etc.. The low 
energy physics of the A^ = phase is identically that of 
an insulating SDW phase with perfect nesting FS. 

The actual field values for which the MR mode drowns 
in the continuum in phase A^ = cannot be taken 
too seriously, since for that range of field values, the 
approximation 2So/hLUc ^ 1 breaks down completely. 
However, there is little doubt that the behaviour found 
above holds qualitatively since the approximation be- 
comes good again at higher fields. 



IV. CONCLUSION 

In addition to proving that within the RPA, more than 
one MR mode is usually alive in the single particle gap, 
at g|| = 1 or 2 (possibly 3, or 5,6,8, etc. due to the 
vicinity of Ribault phases) and computing their ener- 
gies as a function of field, we have derived for the first 
time numerical estimates of the very large effective mass 
anisotropy of the MR. We have described how, as the 
magnetic field increases, UQC phases at small quantum 
numbers progressively lose their MR and, above a thresh- 
hold field, retrieve a low energy excitation spectrum of a 
conventional SDW with perfect nesting. Our results are 
obtained within the RPA, in the "very weak" coupling 
limit {2Sn/^0Jc ^ 1)- In Bechgaard salts, this ratio is of 
order .4 below 6 T, and increases with field. Between 8 
T and 20 T, it is actually larger than 1, then it decreases 
and becomes small at very large fields. Thus, significant 
corrections may be expected, in a weak coupling RPA, 
(i. e. A <C 1) to all estimates derived here, in particular 
between 8 T and 20 T. We have found that the lowest 
MR energy minima of the positive UQC phases could 
be sometimes below Sn. Those values might be signif- 
icantly lowered by self-trapping effects |jl^, which are 
outside the scope of linear response theory |l^. How- 
ever the presence of MR with XN^m < -5 does not seem 
to be garanteed for all A^. The effective masses should 
increase due to MR-MR scattering processes, etc.. The 
importance of those effects will be gauged by comparing 
to the estimates given here. The latter should be a useful 
guideline for experimental efforts to determine the MR 
parameters. 
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FIG. 3. Dispersion of the MR modes in the transverse di- 
rection for the two low energy modes in UQC phase 5; results 
arc similar for UQC 4. The effective mass m± is found to be 
about 10^ larger than m||. 

FIG. 4. MR energy minima as function of H for the high 
field UQC phases = 1 and A'^ = 0. The phase 1 loses its high 
energy MR above H = 7.5T. The phase loses its MR one 
by one as the field increases. Above H ~ 12. 5T, the only col- 
lective modes left are the Goldstone bososns connected with 
the SDW broken symmetries 
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FIG. 1. Network of TN±miq±) lines in the {T,H) plane 
for 6T< H < 12.5 T. Dotted lines are the upper limit of the 
family obtained by varying q±. At T = OK, 26n — 3.53 x Tcn 



FIG. 2. MR energy minima as function of H for the generic 
UQC phase N = 4 and neighbours. The lowest energy MR is 
for =1, the second lowest one for 2. Both are almost H 
independent. Both modes are discontinuous at the transition 
from sub-phase to sub-phase. The MR with g|| =6, due 
to contamination from the Ribault phase —2 lies highest in 
energy, appears only for H > 3.75T and has a monotonic 
negative variation with H. A mode at g|| =3 appears very 
close to the single particle continuum in phase 5 at J? ~ 3.6 
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